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ABSTRACT 

The amplification of a magnetic field due to the Richtmyer-Meshkov instability (RMI) is investigated 
by two-dimensional MHD simulations. Single-mode analysis is adopted to reveal definite relation 
between the nonlinear evolution of RMI and the field enhancement. It is found that an ambient 
magnetic field is stretched by fluid motions associated with the RMI, and the strength is amplified 
significantly by more than two orders of magnitude. The saturation level of the field is determined 
by a balance between the amplified magnetic pressure and the thermal pressure after shock passage. 
This effective amplification can be achieved in a wide range of the conditions for the RMI such as 
the Mach number of an incident shock and the density ratio at a contact discontinuity. The results 
suggest that the RMI could be a robust mechanism of the amplification of interstellar magnetic fields 
and cause the origin of localized strong fields observed at the shock of supernova remnants. 
Subject headings: instabilities — MHD — magnetic fields — shock waves — turbulence 



1. INTRODUCTION 

The Richtmyer-Meshkov instability (RMI) is of cru- 
cial importance in a variety of applications including 
astrophysical phenomena and laborato ry experiments 
(|Brouiliettel 120021 ; iNishihara et~aTll2010D . The RMI oc- 
curs when an incident shock strikes a corrugated contact 
disco ntinuity separating two fluids wi th different densi- 
ties (|Richtmverlll960t iMeshkovi [1969T ) . Because of the 
corrugation of the interface, the surface profiles of the 
transmitted and reflected shocks are also rippled. The 
RMI is driven by the vorticity left by these rippled shocks 
at th e interface and in the fluids ([Wouchuk fc Nishihara! 
[1991 . 

Supernova remnants (SNRs) are expected as a site 
of magnetic field amplification. Thin shell structure 
of the front shock in young S NRs observed by hard 
X-rays (e.g., lLong et all 12001 iBamba et all [20031 ) re- 
quires the existence of strong fields with ~ 100 fiG, 
which is estimated from the syncrotro n cooling length 
(jBerezhko et al.ll2003t fVolk et al.ll2005f ). Recent discov- 
ery of the year-scale variability in the synchrotron X-ray 
emission at the downstream regions of supernova shocks 
suggests that the magnetic field should be amplified up 
to the level of milligauss the re (jUchivama et al.l l2007t 
lUchivama fc Aharonianll2008t) . 

In this paper, we consider the origin of this localized 
milligauss field in SNRs. Since the typical magnetic field 
in the int erstellar medium (ISM) is of the order of mi- 
crogauss (lBeckll2001l : lHeiles fc Trolandl[2005l) . the ampli- 
fication beyond simple shock compression is necessary to 
achieve the milligauss-order magnetic field. It is reported 
that observed synchrotr on images could b e modeled by 
random magnetic fields (jBykov et al.ll20081 ). but the am- 
plification mechanism of the field is st i ll an open ques- 
tion. As discussed in lUchivama et all (|2007t ). the pres- 
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ence of milligauss magnetic field has a critical meaning in 
the long-standing para digm of cosmic-ray proton accel- 
eratio n in young SNRs (|BeUlll978l iBlandford fc Ostrikerl 
Il978f l. The investigation on the physical origin of such 
strong field is then an inevitable step to clarify the role 
of the milligauss field. 

It is widely known that the ISM has nonuniform 
density structures of warm and co ld neutral media 
formed by the thermal i nstability (Field et al.l 119691 : 
iKovama fc Inutsukal [200l llnoue et alll2009f ). The elec- 
tron density fluctuation in hot ionized medium is also 
observed as an eviden ce of the interstellar turbulence 
(|Armstrong et al.lll995f) . The RMI can take place when 
supernova shocks pass through these inhomogeneous ma- 
terials, and which is a quite common event in the ISM. 
The low- and mid-temperature parts of the ISM is par- 
tially ionized by cosmic rays. Although the ionization 
fraction is very small, it is well described by the ideal 
MHD approximation. Hence the interstellar magnetic 
fields could be affected by the fluid motions of the RMI. 

Turbulence and magnetic field amplification has been 
studied for shoc k wave propagation throu g h the inho- 
mogeneous ISM dGiacalone fc Jokipiil l2007t llnoue et"aTI 
120091 120121 ; IGuo et al.l I2012T). Field amplification is in- 
ferred to take place both in the star's collapsing core 
(jEndeve et all 120101) and in the progenitor wind over- 
run by the shock (jChevalierl I1992T ) , and in both cases 
the shock itself has been sugg e sted a s the cause of 
the amplification, llnoue et al.1 (|2012h examined the 
three-dimensional magnctohydrodynamical interaction 
between clumpy interstellar clouds formed by the ther- 
mal instability and a strong shock wave of young SNRs 
with the typical velocity ~ 3000 km s _1 . It is found 
that the maximum strength of the interstellar field in- 
creases rapidly owing to the shear layer at the cloud in- 
terface, in addition to rather gradual growth of the aver- 
age field strength by the so-called turbulent dynamo (or 
turbulent stretching of magnetic field). The maximum 
field strength reaches 1 mG in their simulations, which 
is larger by an order than the volume averaged strength. 
The strong-field regions are localized in thin layers at the 
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interface. 

Turbulent dynamo is one of the promising effects to 
amplify the magnetic field in the postshocked regions. 
This process dominates the growth of the volume aver- 
aged fiel d strength, which has been recognized by man y 
authors (jGiacalone fe Jokipii|[2007t IMizuno et al.ll201lD . 
On the other hand, the size of the strong-field regions 
inferred from X-ray o bservations is typically 0.05 pc 
(|Uchivama et al-lfeOOTt ). and this would be much smaller 
than the preshock inhomogeneity scale ~ 1 pc. Such lo- 
calized strong field should be attributed to the shock- 
cloud interaction, and characterized by the maximum 
strength rather than the average. Therefore, the evo- 
lution of the maximum field strength is a key to under- 
stand the formation of the localized milligauss regions. 
In this paper, we demonstrate the RMI could be an ef- 
ficient mechanism of the field amplification at the cloud 
interface, which can enhance the maximum strength lo- 
cally by a factor of more than 100. 

Shock-cloud interaction with a magnetic field has been 
studi ed numerically by multi-dimensional simulations 
(e.g., Ivan Loo et all 120071 : iShin et all l2008t ). In their 
works, the time evolutions of the entire region of a round- 
shaped cloud are considered. In contrast, we focus on the 
hydrodynamical instability arising from the spacial fluc- 
tuations at the cloud surface aiming to see much smaller 
scale structures. The advantage of this approach is that 
we can make use of the quantitative knowledge of RMI 
for the interpretation of the results. 

MHD effects on the RMI are investigated firstly in 
the context of the dy namics of the m agnctosphcrc 
(|Wu fe Roberts! Il999f ). iSamtanevI (|2003t) has demon- 
strated by numerical simulations that the growth of the 
RMI is suppressed in the presence of a magnetic field. 
Depending on the field direction relative to the incident 
shock surface, the magnetic field works differently on 
the fluid motions of RMI. When the magnetic field is 
aligned with the motion of the shock, the suppression of 
the instability is caused by changes in the shock refrac- 
tion process at the contact discontinuity. The role of the 
magnetic field is to prevent the deposition o f circulation 
on the interface (jWheatlev et al.ll2005L r2009Q . When the 
magnetic field is parallel to the shock s urface, the insta - 
bility is stabilized by the Lorentz force (|Cao et al.ll2008t ). 

Previous works on the MHD RMI is mostly focused on 
the suppression effects by a strong magnetic field. Here, 
we concentrate our discussions on the evolution of an am- 
bient magnetic field affected by the nonlinear motions of 
the RMI. The velocity of the supernova shocks is typi- 
cally 100-1000 km s , and then the Mach number can 
be 10-100 in the warm diffuse gas. The thermal instabil- 
ity produces the density gap between the cold and warm 
neutral media and the ratio is about 10-100. Thus, the 
RMI associated with the supernova shocks is character- 
ized by a large Mach number and a large density ratio. 
In this paper, we show the results of MHD simulations 
performed with a variety of initial conditions including 
such extreme cases. 

Single-mode analysis of the RMI is adopted in order 
to examine the field evolution and understand the phys- 
ical mechanism. The plan of this paper is as follows. In 
§2, the basic equations, initial conditions, and numerical 
method are described. In this section, we make a brief re- 
view on the characteristics of the RMI, which is useful for 



the interpretation of our results. Numerical outcomes of 
two-dimensional MHD simulations are shown in §3. The 
time evolutions of the field for various initial conditions 
are examined in this section. In §4, the saturation mech- 
anism of the amplified field and parameter dependence 
of the results are discussed. Finally, our conclusions are 
summarized in §5. 



2. METHOD 

2.1. Basic Equations and Initial Conditions 

To study the nonlinear evolutions of the RMI, 
equations of ideal MHD are solved; 
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where p, v, and B are the mass density, velocity, and 
magnetic field, respectively, and e is the total energy den- 
sity, 
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In our simulations, the equation of state for ideal gas is 
adopted with the adiabatic index 7 = 5/3. 

The initial configuration for our single-mode analysis 
is illustrated in Figure [TJ1. Two fluids with different den- 
sities, p\ and P2{> Pi), are separated by a contact dis- 
continuity at y = y c d- A shock propagating through 
the light fluid "1" (y > F c d) hits the interface at a time 
t = 0. The incident shock velocity is Ui(< 0) and the 
fluid velocity behind the shock is V\(< 0). Let Pq be 
a uniform pressure of the both fluids before shock pas- 
sage. The Mach number of the incident shock is defined 
as M = \Ui\/c s \ where c s \ = (7P0/ Pi) 1 ^ 2 is the sound 
speed of the fluid "1" . Here the x- and y-axis are set to 
be perpendicular and parallel to the shock surface. 

The interface has an initial corrugation of a sinusoidal 
form, y = Y cc i + iJjq cos(kx), where tpo is a corrugation 
amplitude, k = 2ir/\ is the perturbation wavenumbcr, 
and A is the wavelength. This spatial corrugation of the 
interface is an essential ingredient for the occurrence of 
the RMI. 

The initial geometry of a magnetic field is assumed to 
be uniform with the size of \B\ = Bq in the preshocked 
region (y < Y"j s ). As for the field direction, three cases 
are considered and those are for the cases of perpendic- 
ular MHD shocks (B x = Bq), parallel shocks (B y = Bo), 
and oblique shocks (B x = By). The physical quanti- 
ties in the postshocked region behind the incident shock 
are calculated from the Rankine-Hugoniot conditions for 
MHD shocks. 

The initial configuration depicted by Figure \Va can be 
characterized by only four non-dimensional parameters. 
The sonic Mach number M parameterizes the incident 
shock velocity. The contact discontinuity is expressed 
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by the density jump P2I P\ and the ratio of the corru- 
gation amplitude to the wavelength tpo/X. The initial 
field strength is given by the plasma beta /3o = 87rFo / Bq 
which is the ratio between the gas and magnetic pres- 
sures defined at the preshocked region. Thanks to the 
simplicity of the system, a variety of situations can be 
examined only by choosing these four parameters; M, 
P2/ Pi, 4'o/\ and fa. 

2.2. Characteristic Scales for the RMI 

Figurc[lJ> is a schematic picture around the contact dis- 
continuity after the shock passage. Now, the contact dis- 
continuity moves with the velocity v*. The pressure P* 
becomes higher compared to the initial Pq. The incident 
shock is transmitted into the heavy fluid "2" (y < Y c ^) 
which moves with the velocity Ut(< 0), and the density 
behind the transmitted shock is p\ . Similarly, a reflected 
shock forms with the velocity U r (> 0) and compresses 
the fluid "1" to the density p*. If the interface is cor- 
rugated, the shock fronts of these two waves are also 
rippled. Then, refraction of the fluids across the shocks 
occurs as indicated by the arrows in Figure [lj). 

Consider fluid motions near the shock surface in a 
frame moving with the transmitted or reflected shock. 
For the both shocks, the upstream fluid moves along the 
y-axis. However, the downstream motion should have the 
tangential component as a consequence of the obliqueness 
of the shock front. When the corrugation amplitude is 
small, the velocities of these tangential flows are given 
by 

Svi = kMv* ~Vi), (6) 

Sv% = k?p t v* (7) 

(jWouchuk fc Nishiharal[l997i) . where Vv = ip (l-U r /Ui) 
and ipt = ^0(1 — Ut/Ui) are the initial ripple amplitudes 
for the reflected and transmitted shocks. The tangential 
shear motion across the contact discontinuity generates 
the vorticity at the interface, and which is the driving 
source of unstable growth for the RMI. 

In the RMI, the shock acceleration is impulsive and 
cause s the pertu r batio n amplitude to grow linearly in 
time. iRichtmven (119601 ) has proposed a simple estimate 
of the growth rate of the amplitude; dip/dt = kip^A*v* 
where V>o = ^o(l — v* /Ui) is the interface amplitude just 
after shock passage and A* = (p* — p%)/(pl + p%) is the 
Atwood number of the postshocked interface. This pre- 
scription is obtained by a generalization of the Rayleigh- 
Taylor formula with an impulsive acceleration. 

More complex analytic formul as for the asymptotic 
growth rate of RMI are derived bv lWouchuk fc Nishiharal 
(|1996l IT9971 ). In the weak shock limit, the linear growth 
velocity i>ii n can be written using the tangential velocities 
Svl and 8v% as 

p\5vl - p* 2 b~vl 

Vlin = pT+pI ' (8) 

and then the timescale of the RMI is characterized by 

kvn n 

Within the linear regime, the spike and bubble grow with 
the same velocity vu n . For the weak shock cases, it is 



sufficient that the vorticity only at the interface is con- 
sidered for the derivation of un n . However, for the strong 
shock cases, the bulk vorticity left behind the rippled 
transmitted shock cannot be ignored. It is known that 
the bulk vorticity reduces the growth of the RMI. There- 
fore, in the strong shock limit, the growth velocity given 
by Equation (El) will be overe stimate by a factor of a few 
(|Wouchuk fc Nishiharal IT997D . Furthermore the growth 
rate of the RM I could be reduced in the p resence of a 
magnetic field (fWheatlev et all 120051 [200$ }. However, 
for sake of simplicity, we used Equations flg]) and ((9|) as 
the typical scales of the RMI throughout our analysis in 
this paper. 

Once the initial conditions are set, the growth veloc- 
ity fun can be derived by solving an appropriate MHD 
Riemann problem. The analytic formula given by Equa- 
tion ([5} intimates that the faster growth can be realized 
by the larger amplitude ifto and/or the larger Mach num- 
ber M. When M > 2, the growth velocity is roughly 
proportional to M. This means that the growth time 
of the RMI is quite different depending on the inci- 
dent shock strength. In the limit of high density ratio, 
Pi I Pi S> 1, the growth velocity is approximately esti- 
mated as ~ fc^ |w*| ~ &?/'o(P2/pi)~ 1 / 2 |Vi|. However, 
within the range of pi/ p\ ~ 1.5-100 considered in our 
analysis, the growth time t rra has little dependence on 
the density jump. 

Just for reference, we estimate the characteristic quan- 
tities related to the RMI occurring in the ISM. The shock 
velocity of SNRs is 100-1000 km s" 1 , so that the typ- 
ical Mach number is M ~ 10-100 in the warm neu- 
tral medium. It is a common feature in the ISM that 
nonuniform structures consisted of warm and cold neu- 
tral media. The cold and dense filamentary clumps and 
the surrounding warm diffuse gas are form ed as a nat- 
ural consequence of the th ermal instability (|Fieldl [19651 : 
iKovama fc Inutsukal l2000f ) . The density ratio between 
the two phases is about 10-100. The length of the fil- 
ament (~ 1 pc) is roughly given by the most unsta- 
ble scale of the thermal insta bility, and the thi ckness of 
the filament is about 0.1 pc (|Inoue et al.l I2009T) . At the 
boundary between the two phases, there exists a tran- 
sition layer, the size of which is determined by the so- 
called Field length (~ 0.01-0.1 pc). It is known that 
the trans ition layer is unstable under corrugational defor- 
mations (jlnoue et alll2006HStone fc Zweibell 120091 ) . The 
most unstable wavelength of this instability at the evap- 
oration/condensation front is typically 0.1 pc and the 
growth time is approximately equal to the cooling time. 
Then we can assume the initial corrugation amplitude for 
the RMI is of the order of the ratio of the transition thick- 
ness to the most unstable wavelength, i.e., ipo/\ ~ 0.1. 

Using these typical quantities of the ISM, the growth 
velocity is estimated as i>ii n /c s i ~ 1-10 and the timescale 
of RMI is roughly given by t rra /t s \ ~ 0.1-0.01, where c s \ 
is the sound speed in the diffuse gas and t s \ = A/c s i 
is the sound crossing time. Assuming c s i ~ 10 km s _1 
and A ~ 0.1 pc, these values are corresponding to vy m ~ 
10-100 km s" 1 and t tm - 10 2 -10 3 yr. Then the growth 
time t TTa is much shorter than the cooling time £ nr >oi w hich 
is a few Myr in the typical ISM (|Inoue et all 120091 ). In 
the analysis below, we choose M = 10, pi/ pi = 10, and 
^o/A = 0.1 as the parameters for a fiducial case. 
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2.3. Numerical Scheme and Grid 

We solve the ideal MHD equation s by using a con - 
servative Godunov-type scheme (e.g.. iSano et al.lll998[ ). 
Operator splitting algorithm is adopted in our scheme. 
The hydrodynamical part of the equations is solved by 
a second-order Godunov method, using the exact so- 
lutions of a simplified MHD Ricmann problem. The 
one-dimensional Ricmann solver is simplified by in- 
cluding only the tangential component of a magnetic 
field. The characteristic velocity is then the magne- 
tosonic wave alone, and then the MHD Riemann prob- 
lem can be solved in a way similar to the hydrodynami- 
cal one (jGolella fc Woodward! Il98 l. The piecewise lin- 
ear distributions of flow quantities are calculated with 
a monotonicity constraint following van Leer's method 
(jvan Leer! Il979f h The remaining terms, the induction 
equation and the magnetic tension part of the equa- 
tion of motion, are solved by the consistent MoC-CT 
method (jStone fc Normanl fT99l iClarkd fl996l) . guaran- 
teeing V • B = within round-off error throughout the 
calculation (jEvans fc HawlevllT988h . 

The numerical scheme includes an additional numerical 
diffusion in the direction tangential to th e shock surface 
in or der to care the carbuncle instability ([Hanawa et al.l 
120081 ). The numerical viscosity is added only in the re- 
gions where the characteristics of either fast or slow wave 
converges, i.e., in the regions potentially dangerous to the 
carbuncle instability. This treatment could be crucial in 
our analysis because this numerical instability influences 
the size of baroclinic vorticity, which is generated by the 
misalignment of the pressure gradient of the shock and 
the local density gradient across the interface. 

Calculations are carried out in a frame moving with 
the velocity v* which is the interface velocity after the 
interaction with the incident shock. For convenience, we 
define the locus of y — as where the incident shock 
reaches the contact discontinuity at t = 0. Then the 
postshocked interface stays at y = if it is not corrugated 
initially, or stable for the RMI. The simulations start 
before the incident shock hits the corrugated interface 
(\Yi S — Y" c d| > tpo), and thus the initial time is negative 

t = -IKfc-Kcl/IE/il. . 

The system of equations is normalized by the initial 
density and sound speed of the light fluid "1" , p\ = 1 and 
c s i = 1, and the wavelength of the density fluctuation 
A = 1. The sound crossing time is also unity in this unit, 
t s i = X/c s i = 1. Most of the calculations in this paper 
use a standard resolution of = A y = A/256 unless 
otherwise stated. A periodic boundary condition is used 
in the x-direction, and an outflow boundary condition is 
adopted in the y-direction. The size of the computational 
box in the x-direction is always set to be L x = A. The 
y-lcngth, on the other hand, is taken to be relatively 
larger, L y > 10A. The choice of L y is depending on 
the initial model parameters. All the calculations are 
stopped before the transmitted shock reaches the edge 
of the computational domain. 

3. RESULTS 

3.1. Magnetic Field Amplification 

Single-mode analysis of the RMI and the amplification 
of a magnetic field arc investigated by two-dimensional 
MHD simulations. Figure [2] shows the two-dimensional 



images of the density and field strength at the nonlinear 
regime of RMI. The initial parameters of this fiducial 
model are the Mach number M — 10, the density jump 
Pit 1 Pi = 10, and the corrugation amplitude ipo/\ = 0.1. 
A uniform magnetic field parallel to the shock surface, 
(B x ,B y ) = (Bq,0), is assumed initially, and thus this 
case is for a perpendicular MHD shock. 

In order to elucidate passive evolution of the field, the 
initial strength is supposed to be very weak, /3 = 10 8 
at the upstream of the incident shock. Even though this 
weak field does not affect the dynamics of the RMI, the 
spatial distribution of the field can be modified dramat- 
ically by the fluid motions driven by the RMI. Here- 
after, we use the time normalized by the timescale of 
RMI, t lm — (fcwiin) -1 . The ratio between i rm and the 
sound crossing time t s \ depends on the model parame- 
ters. Most of the cases, i rm is shorter than t s \. For exam- 
ple, t rm = 0.116 for the fiducial model, while t s \ is always 
unity in our simulations. The snapshot data shown by 
Figure [2] are taken at the normalized time kvn n t = 10. 

The density profile shown in Figure [SJi exhibits a 
mushroom-shaped spike as a result of the growth of RMI. 
This figure is a close-up view focused around the contact 
discontinuity at y = 0. The loci of the transmitted and 
reflected shocks at this time are further out of this fig- 
ure, that is, y = — 1.6A and 5.6A, respectively. The spike 
height reaches y ~ 0.7A, which is 7 times larger than 
the initial amplitude i/jq. In Figure [3ji, the time history 
of the distance from the spike top to the bubble bottom 
dsb = \ys — j/t> | is shown as a function of the normalized 
time kv\i n t, where y s and yb are the y-coordinate of the 
spike and bubble. 

Figure [3]j also shows the growth velocities of the RMI, 
which are evaluated by the advection velocities of the 
spike v s and bubble «&. The growth velocity of the spike 
v s increases rapidly just after the incident shock hits the 
interface, and takes the maximum when kvu n t ~ 2. The 
maximum is comparable to the analytic prediction vy m 
given by Equation ([8]). The growth velocity i>i; n is sub- 
sonic compared to the sound speed in the shocked re- 
gions, v lin = 0.20c*! = 0.48c* 2 where c* sl = (7P* / ' p\) 1/2 

and c* 2 = {lP* / P?) 1 ^ 2 ■ but it is supersonic relative to 
the preshocked sound speed, i>i; n = 1.37c s i. The spike 
continues to grow with the velocity about v s ~ 0.2wii n 
until at least kvy m t = 30 in this model. 

The bubble growth, on the other hand, is quenched 
at the earlier phase of the evolution. Th is is a general 
feature of the RMI (jMatsuoka et al.ll2003l ) , especially for 
the strong shock cases as examined in this paper. The 
suppression is caused by the bulk vorticity comes from 
the deformed transmitted shock. It is worth noticing 
that weak spike-like patterns can be seen in the density 
distribution behind the transmitted shock in the heavy 
fluid "2" . This is also associated with the bulk vortic- 
ity generated by the propa gation of a rippled shock (e.g., 
Ilshizaki fc Nish ihara 19971). As the shock separates away 
from the interface, their corrugation amplitude oscillates 
and decreases with time. The amplitude of the bulk vor- 
ticity is larger for the cases with higher shock strength, 
and it could affect the nonl i near behavior of the RMI 
(|Wouchuk fc Nishiharal[l997tlNishihara et al.ll2010Tl . Al- 
though the effects of the bulk vorticity might be impor- 
tant for the interstellar strong shocks, its quantitative 
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details is beyond the scope of this paper. 

Figure depicts the spatial distribution of the field 
strength. The magnetic field is already amplified more 
than 200 times at kv\i n t = fO. The strong field regions 
are localized at the mushroom cap and form thin filamen- 
tary structures along the interface. Figure [3l> shows the 
time evolution of the maximum field strength |-B| max in 
the computational domain. Before the incident shock 
hits the corrugated interface, the y-component of the 
field is nothing in this model. The Mach number M = 10 
is large enough that the density behind the incident shock 
is enhanced by a factor of (7 + l)/(7 — 1) ~ 4 due to 
shock compression. Then the initial |-B x |max is larger by 
the same factor than Bq. 

During the interaction between the incident shock and 
interface, B y appears near the rippled interface due to 
the refraction of the flow. The maximum value l-B^lmax 
grows linearly in time at the very beginning. Then 
each component of the field increases exponentially until 
kv\i n t ~ 2 during the growth velocity of the spike v s is 
comparable to vu n (see Figure[3ji). At this stage, |-Ba:|max 
is still dominant over l-Bylmax- But the both components 
become comparable and evolve in a similar manner at 
the later evolutionary stage kvn n t > 5. The early growth 
rates ob obtained by using a fitting formula, 



(t) oc exp 



( 

\ ^rm 



(10) 



are ob = 1-0, which is shown by the dotted line in the 
figure. The maximum strength increases with time even 
at kvu n t = 30 and exceeds |-B| max ~ 10 3 Bq at the end. 
These results are obviously suggesting that the RMI is a 
quite efficient mechanism of magnetic field amplification. 

3.2. Physical Mechanism of the Amplification 

In this subsection, we consider the physical mechanism 
of the field amplification associated with the RMI. The 
induction equation for the ideal MHD can be rewritten 

as 



f = -(,-V)B 



(B ■ V)v - B(V ■ v) , (11) 



where each term of the right-hand-side stands for advec- 
tion, stretching, and compression. Among these terms, 
net increase of the field can be done only through the 
stretching and compression. 

Figure IU2 shows the magnetic field lines in the neigh- 
borhood of the interface for the fiducial model at kvn n t = 
2. The gray color denotes the higher density part in 
the fluid "2" bounded between the contact discontinuity 
and transmitted shock front. As mentioned above, the 
magnetic field is amplified efficiently at the interface be- 
tween two fluids. The positions of the largest strength 
are shown by the crosses in this figure. The fluid motion 
excited by the RMI stretches the surface area of the con- 
tact discontinuity. It is found that the field amplification 
is predominantly originated from this stretching effect. 

The relative importance of stretching and compression 
can be evaluated from our numerical results. Multiply- 
ing the magnetic field vector B to Equation (fTTj). the 
following scalar equation can be derived; 

\^- t \B\ 2 = -B-(vV)B+B-(B-V)v-\B\ 2 V-v , (12) 



where the last two terms represent stretching and com- 
pression. 

Using the snapshot data same as in Figure^, we cal- 
culate the sizes of the stretching and compression terms 
in Equation (|12j) . which are shown in Figures © and 
|4k, respectively. The same color contours are used in 
these two figures. Focusing on the mushroom-cap region, 
the stretching term takes the maximum, while the com- 
pression term is almost nothing there. Because the fluid 
motions in the RMI is mostly incompressible, so that 
the contribution of compression toward the field ampli- 
fication is negligibly small compared to that of stretch- 
ing. The compression term is larger at the transmitted 
shock front, but that is simply because of shock com- 
pression. The strong field regions along the interface 
perfectly match where the stretching term is dominant 
over the other terms. This is a clear evidence that the 
stretching effect is the major source of the field amplifi- 
cation. 

To verify this interpretation, we compare quantita- 
tively the stretching rate of the interface and the growth 
rate of the field. The stretching rate can be ob- 
tained numerically fro m a nonlinear vortex sheet model 
(jMatsuoka et al.l [20031 ). After the transmitted and re- 
flected shocks have traveled a distance larger than the 
fluctuation wavelength A, the system can be regarded 
as incompressible and irrotational except for the inter- 
face on which nonuniform vorticity is induced by the 
shocks. We treat the interface as a curve with the use of 
a Lagrangian parameter in the x-y plane. The govern- 
ing equations for the nonlinear analysis are the Bernoulli 
equation and kinematic boundary conditions. 

The Bernoulli equation, i.e., the pressure continuous 
condition at the interface, is given by 



Pi 



1 



at 



(V0i 



l>2 



1 



dt 



(V0 2 ) 2 



(13) 



where cf>i (i — 1,2) is the velocity potential defined as 
= Vi and Vi is the velocity of the fluid "i" . In order 
to calculate the interfacial motion as a vortex sheet, we 
need to rewrite Equation (|13[) into an evolution equation 
for the vortex sheet strength k. Here, k = \n\ is defined 
by the circulation r = 4>\ — <f>2 as n = Vr. The detailed 
derivation of the evolution equation associated with the 
vortex-induced velocity is explained in the appendix. 

The obtained shape of the interface at kvn n t = 2 is 
shown by Figure [4]1 The line color indicates the stretch- 
ing rate <7i n t at each Lagrangian point on the interface, 
which is defined as 



Cint 



trm dl 

~T~dl 



(14) 



where I is the line element along the interface. The vortex 
sheet result reproduces surprisingly well the spatial dis- 
tribution of the stretching size in our simulation shown 
by Figure [4j>. The large stretching rate can be seen at 
the spike top, and this feature coincides exactly with 
the strong field region as shown in Figure 0h- The vor- 
tex sheet model predicts that the largest stretching rate 
appears at the side of the mushroom cap in the later 
phase. This is also consistent with the MHD results at 
kvn n t =10 depicted by Figure [2j>. Furthermore the max- 
imum stretching rate <7i n t ~ 0.8 is fairly close to the 
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growth rate of the magnetic field as ~ 1.0. This is an- 
other fact to convince us of the strong relation between 
the field amplification and the stretching effect. 

3.3. Dependence on the Initial Field Geometry 

The orientation of interstellar magnetic fields has 
no correlation with the direction of supernova shocks. 
Therefore, it is interesting to examine the effects of 
the initial field direction on the amplification process. 
A parallel shock case with the initial field (B x ,B y ) = 
(0, -Bq) and an oblique shock case with (B x ,B y ) — 

(Bq/VZ, B0/V2) are performed for this purpose. The 
other parameters are identical to the fiducial model, so 
that the magnetic field is again assumed to be very weak, 
/3o = 10 8 . The physical quantities in the downstream of 
the incident shock are calculated from the jump condi- 
tions for MHD shocks, and the fast shock condition is 
used for the oblique shock case. 

The field amplification due to the RMI is found to be 
independent of the direction of the ambient field. Fig- 
ure O shows the time evolutions of the maximum field 
strength for the models with different initial field ge- 
ometries. First, let us take a look at the early stage 
of the evolution until kvu n t w 10. For all the models, the 
magnetic field increases exponentially at this stage and 
is amplified up to about 100 times as large as the ini- 
tial strength Bo. The stretching term always contributes 
the field amplification dominantly. The difference in the 
maximum strength at t = is caused by the compression 
of the tangential field B x behind the incident shock. For 
the strong shock limit, the ratio of the initial |-B K | max in 
our models should be (parallel) : (oblique) : (perpen- 
dicular) = 1 : 2V2 ■■ 4. The highest field at kvn n t = 10 
is achieved in the perpendicular shock case, but this is 
mainly because of the difference in the initial |_B| max . 
The amplification factor from the initial maximum field 
is rather similar for all these models. The reason is that 
the evolution of the RMI is unaffected by the weak field, 
and then the stretching motion is almost identical no 
matter which direction the ambient field takes. 

Figures [6jx and \6$> show the spatial profile of the field 
strength and the field lines, respectively, for the oblique 
shock case. These snapshots are taken at kvy m t = 10 
the same as Figure [2j The density distribution is al- 
most identical to that in the fiducial model shown by 
Figure |2Ji. Although asymmetric structures can be seen 
in the amplified field because of the shock obliqueness, 
the strong-field regions are concentrated along the in- 
terface at the mushroom-shaped spike. The maximum 
field appears near the top of the spike and the strength 
is I -B I max /-Bo = 192 at this time. This result strongly 
supports that the field amplification in the oblique shock 
case is also caused by the surface stretching of the con- 
tact discontinuity as well as in the perpendicular shock 
case (fiducial model). 

At the later evolutionary stage 10 < kv\i n t < 30, the 
maximum \B\ keeps increasing with time in the perpen- 
dicular and oblique shock cases (see Fig. EJ. When the 
fluid "1" has a tangential field component, B x , the total 
magnetic flux swept by the spike increases as long as the 
growth of RMI continues. On the other hand, if there 
is only a normal component, B y , the magnetic flux in- 
volved by the RMI is fixed by the initial setting. Then 



the size of the maximum field have to be saturated when 
all the flux is confined in a sufficiently small region. This 
is why the parallel shock case shows the saturation of 
the growth, a flat plateau in the time-profile of |J3| max . 
This interpretation implies the resolution effect on the 
amplification factor. Actually the resolution dependence 
can be seen in some models with a weak initial field, and 
this is discussed later in the next section. 

The time evolutions of the average field strength |-B| avc 
are also shown in Figure [5J Here the average is taken over 
the regions that the pressure is larger than 10Po in order 
to eliminate the preshocked regions. These time profiles 
of the maximum and average field are qualitatively quite 
similar to those in the reali s tic in terstellar simulations 
[e.g., Fig. 5 in llnoue et al.l (|2012t) ]. This resemblance 
reminds that the RMI would be a fundamental process 
of the field amplification in the ISM. However, our re- 
sults are still restricted within the two-dimensional pic- 
ture. Obviously three-dimensional consideration will be 
an important next step for further understandings. 

We found some curious features in the field evolution 
for the parallel shock case. Figures[7]i and© show snap- 
shots of the field lines and the size of the stretching term 
at kvn n t = 6. For this case, the position of the maxi- 
mum field is quite different from that in the other cases. 
The field lines are concentrated at the stem of the mush- 
room rather than the interface or the mushroom cap. 
The heavy fluid "2" near the interface moves horizon- 
tally toward the y-axis just after the shock passage. The 
magnetic field is frozen into the fluid and thus it is gath- 
ered near the y-axis by this motion. The converging flow 
can escape toward the mushroom cap along the field lines 
so that it can be realized as incompressible motion. The 
largely amplified region coincides to where the stretching 
term, B y (dv y /dy) w —B y (dv x /dx), is the highest. In- 
terestingly, the magnetic field is amplified selectively in 
the fluid "2" for the parallel shock cases. For the per- 
pendicular shock cases, the velocity pattern is the same, 
but the converging flow basically moves along the field 
lines and thus no field amplification at the stem part. 

4. DISCUSSIONS 
4.1. Parameter Dependence 

In this section, we discuss the dependence of the field 
amplification on the model parameters such as the Mach 
number M, and the density jump P2/P1, and the initial 
field strength (3q. It is found that the amplification of a 
magnetic field by the RMI occurs in a variety of cases 
with a wide range of the initial parameters. The ampli- 
fication factor is more than 100 for most of the cases. 
Therefore, this process should be a quite common phe- 
nomena, especially in astrophysical shock events with a 
high Mach number. 

Figure [5] shows the evolutions of the maximum field 
strength for various models. The dependence on the 
shock strength is shown in Figure ISJz. The Mach num- 
ber ranging from M = 1.5 to 100 are examined. The 
model parameters except for M are identical to the fidu- 
cial model and the initial field orientation is in the x- 
direction. When M > 3, the amplification factor ex- 
ceeds 10 3 in kv\i n t = 30 as shown in the figure. The 
time history of |B| max in terms of the normalized time 
kvn n t looks quite similar for all the models. The am- 
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plification process seems to be independent of the Mach 
number for this parameter range. This is because the 
stretching rate of the interface in the normalized unit 
has little dependence on the Mach number. Thus, even 
for the weak shock case with M = 1.5, the magnetic field 
can be amplified by more than 100 times in kvn n t ~ 10. 
These results are for the cases of the perpendicular shock, 
but the same conclusion can be obtained for the parallel 
shock cases. Notice that the growth velocity is nearly 
proportional to the Mach number. The normalization of 
the time is £ rm = 1.17t s i for the model with M = 1.5 and 
t tm = 1.15 x 10~ 2 i s i for M = 100. The actual timescale 
of the RMI growth is largely different depending on the 
shock strength. 

Figure [8}) demonstrates the dependence on the density 
jump at the interface. The models of the different density 
ratio, pil p\ = 100, 10, 3, and 1.5 are shown in this figure. 
The other parameters and initial field geometry are the 
same as in the fiducial model. Again, the magnetic field 
is enhanced by more than two orders of magnitude in 
kvy ln t ~ 10 for all the models. 

The amplification factor has a positive correlation with 
pi/ pi- Since the growth velocity Uiin has little de- 
pendence on p2j 1 pi, the growth timescale is compara- 
ble for these cases; t Tm — 0.40i s i for p^j ' pi = 1.5 and 
tim = 0.19i s i for p^l Pi = 100. However, the nonlin- 
ear behavior of the spike is found to be quite different 
at the later stage. The spike height reaches y s « 2.6A 
at kvnnt = 30 for the model with pi I pi = 100, while it 
is only y s ps 0.8A for P2/P1 = 1.5. When the density 
ratio pi I pi is larger, the total magnetic flux swept by 
the spike is larger. This magnetic flux near the interface 
becomes a seed field to be stretched by the RMI. Then 
the maximum strength could be much higher as the ratio 
Pi I Pi is larger. For the parallel shock cases, on the other 
hand, this trend can not be seen at all and the saturated 
amplitude is independent of pi/pi- 

The time evolutions of the average field strength | B | avo 
are also shown in Figure [8] The average is taken over the 
postshocked regions swept by the incident and transmit- 
ted shocks. For all the models, the maximum strength is 
much larger than the average. The strong field regions 
are highly localized along the stretched interface as a re- 
sult of the RMI. 

We performed some models that the ratio pi / pi is less 
than unity. Although the phase of the RMI becomes re- 
versed in those models, the magnetic field is amplified 
by many orders of magnitude in the same manner. It 
should be noted that the growth velocity given by Equa- 
tion © is applicable for the reflected rarefaction wave. 
The phase inversion is expressed by the negative value 
of the growth velocity. The spike can penetrate deeply 
into the light fluid even for the rarefaction cases, and its 
structure does not change by much compared with that 
in the reflection shock cases. The efficient field amplifi- 
cation occurs mostly near the interface of the spike for 
the both cases. 

4.2. Saturation Level of Magnetic Field 

The dependence of the field amplification on the ini- 
tial field strength is shown by Figure |H1 The maximum 
strength |B| max in this figure is normalized by (87TP*) 1 / 2 
where P* is the postshock pressure. The initial plasma 
beta of these models are /3 = 10 8 , 10 4 , 100, and 1. 



The parameters except for /3q are identical to the fiducial 
model and the magnetic field is initially aligned along the 
x-direction for all the models. When /3q is large enough, 
the initial field is amplified more than two orders of mag- 
nitude. However, if /3q < 100, the amplification is lim- 
ited; the factor is ~ 10 for (3q = 100 and only a few for 
A> = 1. 

When the initial strength of the ambient field becomes 
larger, the saturation level at the nonlinear regime is 
almost converged independent of /3q. The upper limit 
of the field strength is about ft* = 8nP* /|-B| 2 lax ~ 10, 
which is roughly equal to the equipartition value to the 
thermal pressure after the shock heating. More precisely, 
the saturated magnetic field strength is determined by a 
balance with the kinetic energy of the growth velocity of 
the spike, l-B^ax/^ 71 " ~ P-ii^Un/^- F° r the strong shock 
cases, both of them can be comparable to the thermal 
energy As seen from Figure Wt>, the amplified field dis- 
tributes along the interface at the mushroom cap with 
the thickness of ~ 0.01A-0.02A. If the magnetic pressure 
becomes comparable to the gas pressure, the interface 
stretching is reduced by the Lorentz force, and the in- 
crease of the maximum field strength will be saturated. 
Then the potential limit of the field amplification is de- 
termined by a condition that the RMI is suppressed by 
the amplified field itself. 

In this figure, higher resolution results arc also shown 
by the dotted curves. The same colors denote the same 
initial conditions. The resolution dependence can be 
seen for the models with a weaker initial field, where 
the higher resolution gives the higher saturation level. 
However, the upper limit of the saturated field has little 
dependence on the grid size. The parallel shock models 
shows qualitatively the same results for the dependence 
on the initial /3q and the grid resolution. Therefore, we 
can conclude that the existence of the upper limit for the 
saturated field is a robust nature. 

Even when the initial field strength is strong ((3q ~ 1), 
the magnetic field can be amplified by a huge factor if 
the Mach number is much larger than about 10. For the 
strong shock cases, the postshock pressure P* is larger 
than Po by many orders of magnitude. Then the mag- 
netic pressure of the initial field, Bq/8it, can be regarded 
as rather weak in terms of P* . Therefore, the vorticity is 
deposited very close to the interface, and then the RMI 
can grow even in the presence of the strong field. As a 
result of the nonlinear growth of the RMI, the magnetic 
field is amplified up to the thermal value of the postshock 
pressure P* . 

The plasma beta in the ISM is u sually close to unity 
HBeckl 120011: iHeiles fc Trolandl 120051) . When /3 - 1 ini- 
tially, the upper limit of the amplification factor is of 
the order of \B\ max /B ~ (P P*/P ) 1/2 ~ (P*/P ) 1/2 - 
If P* /Po > 10 4 , 100-fold enhancement of the magnetic 
field can be expected. The postshock pressure P* is de- 
termined just by the Mach number M and density ratio 
Pi/ Pi- The pressure ratio P* / Pq is nearly proportional 
to the square of the Mach number, but has little de- 
pendence on the density ratio. When the Mach number 
is larger than M > M cr it « 50, the shock heating en- 
hances the pressure by more than four orders of mag- 
nitude P*/Pq > 10 4 . The Mach number of the shocks 
associated with SNRs is sufficiently larger than this crit- 



8 



ical value. If the initial field strength is assumed to be 
about the average value of the ISM, ~ 5 fiG, then the en- 
hanced maximum field reaches to the level of milligauss. 
Therefore, this mechanism could explain the origin of 
the strong field spots observed i n the downstream re- 
gions of the front shock of SN Rs (jUchivama etaLll2007t 
lUchivama k Aharonianll2008[) . 

The critical Mach number M cr ; t « 50 is independent of 
the initial corrugation amplitude of the interface, but the 
size of ipo affects the growth velocity and timcscale of the 
RMI. Assuming ipo/\ = 0.1, the growth time i rm is much 
shorter than the sound crossing time, t rm < 0.03t s i, 
when M > 50. Since the cooling time is comparable 
to t s x, the magnetic field can be amplified significantly 
before the decrease of the pressure by the radiative cool- 
ing. The ionization fraction of the postshocked media is 
sufficiently large enough that the non-ideal MHD effects, 
such as the ambipolar diffusion and ohmic dissipation, 
are ignorable. Then the dissipation timescale of the am- 
plified field is much longer than the RMI timescale. 

In our numerical analysis, a discontinuous density 
jump is assumed at the interface. However, for the cases 
of the ISM, there must be a transition layer with a finite 
thickness between two different media. Then we also 
examined the effects of the transition layer on the field 
amplification. The density profile of the layer is approx- 
imately given by a function cx tanh[(y — Y c d)/L t ] where 
L t denotes the thickness of transition layers. Although 
the growth rate of the magnetic field decreases with the 
increase of the thickness L t , the saturation level of the 
maximum strength is found to be independent of L t if 
the transition layer is thinner than the corrugation wave- 
length L t < A. Thus the field amplification by the RMI 
could be realized even in the typical situations of the 
ISM. 

4.3. Observed Emissions and Particle Acceleration 

If the entire region of SNR RX J1713.7-3946 is mag- 
netized up to B ~ 1 mG, the origin o f gamma-rays 
cannot be leptonic ([Uchivama e t al. 2007]), which is in- 
ferred from the observed power ratio of the synchrotro n 
to inverse-Compton emissions (jAharonian et al.1 12007D . 
However, when the strongly magnetized regions is caused 
by the RMI after the shock-cloud interaction, the aver- 
aged field strength of the overall SNR would be much 
smaller than milligauss, depending on the volume fill- 
ing factor of the cloud in the upstream ISM. Therefore, 
the origin of gamma-ray from RX J1713. 7-3946 should 
be discussed carefully taking account of other observa- 
tional fea tures, such as the spectrum of the gamma-ray 
emission (|Abdo et al .l l2011l; Ilnoue et al.ll2012h and struc- 
ture o f surrounding ISM (jEllison et al.H2010HFukui et all 

120H . 

As we estimated in §2.2, the timescale of field ampli- 
fication by the RMI in young SNRs is about 100 yr, 
which indicates that the high-energy electrons acceler- 
ated at the forward shock (or the transmitted shock) 
by the diffusive shock acceleration are cooled down due 
to the synchrotron loss before the field strength reaches 
milligauss-order. Thus, in order to explain observed 
synchrotron X-ray hot spots, an additional accelerator 
other than the forward shock is necessary to produce 
brightening of the synchrotron X-rays at the region of 
B ~ 1 mG. Ilnoue et"aLl (|2012[ ) discussed that the re- 



flected shock waves induced by the shock-cloud interac- 
tions can accelerate the electrons even in the downstream 
region of the forward shock, and which could explain the 
synchrotron emissions from strongly magnetized regions. 

5. SUMMARY 

We have investigated the evolution of a magnetic field 
associated with the RMI by using two-dimensional MHD 
simulations. In terms of the field amplification, the im- 
portance of "laminar stretching" driven by the RMI at 
the interface is successfully demonstrated. In our single- 
mode analysis, an incident shock propagating through a 
light fluid is considered to encounter a contact surface of 
a heavy fluid. When the interface is spatially corrugated, 
the RMI takes place and a mushroom-shaped structure 
develops in the density profile. An ambient magnetic 
field is initially supposed to be uniform and subthermal. 
Our numerical results for various situations suggest that 
the RMI is an efficient mechanism of the amplification 
of the interstellar magnetic fields. The main conclusions 
are summarized below. 

1. The fluid motions associated with the RMI 
strengthen an ambient magnetic field by many or- 
ders of magnitude. This phenomenon can be seen 
in a wide range of the initial parameters. The am- 
plification factor is almost independent from the 
Mach number of the incident shock and the ini- 
tial field direction, so that it could occur even for 
the cases with a weak shock and/or a small density 
jump. Therefore we can conclude that the RMI is 
a robust mechanism of the ambient field amplifica- 
tion. 

2. The physical mechanism of the field amplification is 
stretching associated with the nonlinear evolution 
of the RMI. The magnetic field is amplified effi- 
ciently at where the stretching term of the induc- 
tion equation is predominant over the other terms. 
In most cases, the strong field regions are localized 
along the mushroom-shaped interface and form fil- 
amentary structures. Curiously, only for the paral- 
lel shock cases, the maximum field appears at the 
stem of the mushroom in the heavier fluid. 

3. The amplified magnetic field is saturated when the 
magnetic pressure becomes comparable to the ther- 
mal pressure after the shock heating. This is be- 
cause of the suppression of the RMI through the 
Lorentz force of the amplified magnetic field. If the 
Mach number of the incident shock is larger than 
about 50, we can expect at least more than 100-fold 
enhancement of the initial field. Thus the RMI can 
be a promising origin of the interst ellar strong fields 
obser ved at the shock of SNRs (jUchivama et al.l 
l2007h . 

The MHD RMI would play an importa nt role not 
only in other astrophysical phenomena (e.g.. Ilnoue et all 
12011 but also in many scientific fields such as interplan- 
etary shocks {Wu 2003 ) and inertial c onfinement fusion 
(jLindl et all 11993 iHolmes et al.lH999l ). In this paper, 
the evolutions of an external magnetic field are examined. 
However, particularly for the case of laser plasmas, a self- 
generated field often cannot be ignorable. The baroclinic 
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term generates a magnetic field as well as the vorticity. 
Then proper treatment of two-fluid effects should be in- 
cluded in the analysis for that purpose, and which will 
be an interesting subject for our future work. 

The RMI has been st udied extensively by labo- 
ratory experiments (e.g. , iNiederhaus fc Jacobs! 120031 : 
iChapman fc Jacobs] 120061) . Laser plasmas can be a 
new platform to examine the RMI in laboratories 
(|Dimonte fc Remington! fl993() and the inclusion of the 



effects of a magnetic field will be possible in such exper- 
iments (|Kuramitsu et al1l2011l ). Therefore the magnetic 
field amplification proposed in this paper could be tested 
by laser experiments in the near future. 

We thank Shu-ichiro Inutsuka, Tomoyuki Hanawa, and 
Takahiro Kudoh for useful discussions and comments. 
Computations were carried out on SX-8R at the Cy- 
bermedia Center and SX-9/B at the Institute of Laser 
Engineering of Osaka University. 



APPENDIX 

VORTEX SHEET MODEL 

Using the circulation L = <fii — <p2 and average velocity potential $ defined by $ = (</>i + </>2)/2, we rewrite the 
Bernoulli equation (|13[) as 

= 2A Aqq-\ k ■ k - aA k ■ q , (Al) 

Dt Dt H H A y v ' 

in which the derivative D/ Dt is given by 

D d an 
_ = V , V = q+ — , 

where q = V<&, k = Vr, A* is the Atwood nu mber defined by A* = (p* - p %)/(p* + 0%), and a = a(A*) (\a\ < 1) is 
a weighting factor such that a ^ for A* =^ dMatsuoka fc Ni shihara 2006). Here, the spatial derivative is taken at 
the interface. 

We regard the interface in the RMI as a curve in the x-y plane, and parameterize it using a Lagrangian parameter 
9. The velocity of the interface (x,y) = (X(9,t),Y(9,t)) is derived as 

X t = U+^ K , Y t =V+^n (A2) 
2s g 2sg 

(|Baker et al.| [l982: Matsuoka & Nishihara 2006|), where the vortex sheet strength k is defined by k = dT/ds = Tg/sg, s 
is arc length of the sheet, and the subscript denotes the differentiation with respect to the variable. The vortex-induced 
velocity U = U(9,t) and V — V(8,t) are given by the Birkhoff-Rott equation: 

TJI0 A = - J- f smh[Y(e,t)-Y(e',t)] K (d',t)sg(9')de' 
[ ' ' At: J_ n cosh [Y(6, t) - Y(8', t)} - cos [X(9, t) - X(9', t)} + S 2 ' { ' 

V(f) t) = ±- r sm[X(e,t)-X(6',t)] K (6',t)s e (e')d0' 
[,> Air J^cosh[Y(6 1 t)-Y(d' 1 t)}-cos[X{9,t)-X(6>,t)} + 6 2 1 ' 

(Birkhoff 1962]; iRottl [19561 ) . whe re we regularize the Cauchy integral using Krasny's S (|Krasnvill987| ). 

Differentiating Equation (|A1|) with respect to 9, we obtain the following Fredholm integral equation of the second 
kind: 

2A* , x (1 + &A*)k, x A*+a, 9 , . , . 

K t - (XgU t + YgV t ) - ± ^^( Xe Ug + ygVg) + — ( K 2 )g . (A5) 

sg Sg Asg 

Solving Equations (| A2|) and (|A5[) simultaneously, we can determine the motion of a vortex sheet in the RMI. 

For the fiducial model, the fluctuation amplitude after shock passed is approximately calculated as V'oM = 0.06188. 
The density for the both layers are p*/pi = 6.716 and p\l p\ = 39.549, so that the Atwood number is A* = —0.7097 
for this case. The initial condition in Figure HJi is then given by 

2 sin 

X(6,0) = 6, Y(6,0) = a cos6 , k(8,0) = Vr" > 

se(0) 

where the initial amplitude ao is set to eto = 2itiPq/\ = 0.3888. 

In the numerical calculation, the stretching rate, i.e., the rate of the temporal change of length s(9) at the Lagrangian 
point 9 is defined by 

1 ds(9) 1 s(0,t + At)-s{6,t) 
s(6) dt s(9) At 1 ' 

where s = [s(t + At) + s(t)] /2 and the length s(8) is calculated by J sgd9, in which the integral is performed as the 
spectral integration with respect to 9. Here, we set the time step At = 0.0002, the regularized parameter 5 — 0.15, 
and the weighting factor a = —0.05. We adopt the trapezoidal rule and the forth-ordcr Runge-Kutta method for the 
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spatial and temporal integration, respectively. For detailed numerical schemes, refer to iMatsuoka fc Nishihara] ([2006D 
and references therein. 
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Fig. 1. — (a) Schematic picture of the initial configuration for single-mode analysis of the RMI. Two fluids are divided by a contact 
discontinuity (CD). The densities of the light fluid "1" and heavy fluid "2" are pi and p2, and the uniform pressure for the both fluids is 
Po- The interface is corrugated sinusoidally with the wavelength A and the amplitude ipo- An incident shock (IS) propagates in the light 
fluid "1" with the shock velocity Ui. The shock strikes the corrugated interface at a time t = 0. Here Vi is the flow velocity behind the 
incident shock, (b) Sketch of the shock-front shapes after the incident shock hits the corrugated interface. The transmitted shock (TS) 
and reflected shock (RS) travel from the contact discontinuity in the opposite direction with the velocities Ut and U r , respectively. The 
pressure and velocity at the contact discontinuity are P* and v* , and the densities behind the transmitted and reflected shocks are p\ 
and pj. Because of the obliqueness of the shock surface, tangential flows, Sv* and Sv^, are generated at the both side of the interface. 
Refraction of the fluid motions at the transmitted and reflected shocks arc shown by the thick arrows. 
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Fig. 2. — Spatial distributions of (a) the density and (b) the magnetic field strength at a nonlinear stage of the RMI for the fiducial 
model. The model parameters are M = 10, pi/ pi = 10, ipo/\ = 0.1, and /3o = 10 8 . The direction of the initial ambient field is in the 
x-dircction, or perpendicular to the incident shock velocity (perpendicular MHD shock). These snapshots arc taken at the normalized time 
k^lin^ — 10- The maximum field strength at this time IS 1 1> | max /S = 272. 
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Fig. 3. — (a) Evolution of the growth velocities of the spike v a and the bubble Vb, and the distance from the spike top to the bubble 
bottom d 8 }, in the fiducial model, which are normalized by the asymptotic linear growth velocity un n and the wavelength of the fluctuation 
A. (b) Time profile of the maximum of the magnetic field strength |B| max in the fiducial model. The maximum values of each component, 
-EMmax and |B a |max, are also shown. The exponential growth at the early stage until kvn n t = 2 can be fitted by a function oc cxp(cr s t/t rm ) 
with (7b = 1-0 where f rm = (kvn n )~ 1 is the characteristic timescale of the RMI. 



13 




Fig. 4. — (a) Magnetic field lines for the fiducial model at kv\j n t = 2. The gray color denotes the higher density regions in the fluid 
"2" compressed by the transmitted shock. The positions of th e m aximum field arc shown by the crosses. (b,c) Relative importance of (b) 
"stretching" and (c) "compression" in the induction equation 1 1 -'[) calculated from snapshot data of the fiducial model at kvy lrL t = 2. Each 
term is normalized by a constant Bokvi in and B = B/\B\ is a unit vector. The same color bars are used for these figures, (d) Interface 
profile predicted by a vortex sheet model. The model parameters corresponding to the fiducial model are used for the numerical calculation. 
The line color indicates the stretching rate of the interface cr- lnt , i.e., the rate of the temporal change of length at each Lagrangian point. 
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Fig. 5. — Time evolution of the maximum field strength |B| max shown as a function of the normalized time kvn n t for the cases with 
different orientations of the initial field. In the perpendicular shock case (fiducial model), the initial field is (B x ,B y ) = (Bq,0). For 
the parallel shock case and the oblique shock case, the uniform field direction is assumed to be (B x ,B y ) = (0, -Bo) and (B x ,B y ) = 
(Bq/\/2, B0/V2), respectively. The other parameters are the same as in the fiducial model. The thin dashed curves are time evolution of 
the average field strength |B| a vc for each case. The average is taken over only the postshocked regions. 
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Fig. 6. — Spatial distributions of (a) the magnetic field strength and (b) the field lines for the oblique shock case. The direction of the 
initial ambient field is 45-degree to the i-axis. The other parameters are identical to those in the fiducial model. These snapshots are taken 
at the normalized time kvn n t = 10. The higher density regions bounded between the contact discontinuity and the transmitted shock front 
arc depicted by the gray color in (b). 
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Fig. 7. — (a) Magnetic field lines of the parallel MHD shock case at kvn n t = 6. The model parameters are identical to the fiducial model, 
but the initial field direction is parallel to the incident shock velocity for this case. The gray color denotes the higher density regions in 
the fluid "2" compressed by the transmitted shock. The position of the maximum field is shown by the cross mark, (b) Distribution of the 
size of "stretching" term in the induction equation 1 11 211 calculated from the snapshot data same as in Figure [TJi. The stretching term is 
normalized by the same constant Bokv\[ n as in Figure [3J>. 
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Fig. 8. — (a) Dependence of the maximum field strength on the Mach number M. The time profiles of |B| max are shown as a function 
of the normalized time kvn n t for the models with different shock strength; M = 100, M = 10 (fiducial model), M = 3, and M = 1.5. The 
model parameters except for M are identical to those of the fiducial model. The thin curves are time evolution of the average field strength 
-B|avc for each case. The average is taken over only the postshocked regions, (b) Dependence of the maximum field strength on the density 
jump pil P\- The time evolution of |B| max are drawn for the models with different density ratio; P2/P1 = 100, P2/P1 = 10 (fiducial model), 
P2/P1 = 3, and pil p\ = 1.5. Other than P2/P1, the same parameters as in the fiducial model are used. The thin curves are time evolution 
of the average field strength |B| a vo for each case. 
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Fig. 9. — Dependence of the saturation level of the amplified magnetic field on the initial plasma beta /3q = 8ttPq/Bq. The maximum 
field strength |B| max is normalized by (SirP*) 1 / 2 using the postshock pressure P* . The model parameters are the same as in the fiducial 
model except for the initial field strength By, which is indicated by the line colors; /3q = 10 8 (fiducial model). 0q = 10 4 , 0o = 10 2 , and 
00 = 1. The initial field direction is assumed to be in the ^-direction for all the models. For the purpose of comparison, higher resolution 



results are also shown in this figure by the thin dashed curves. The grid size of the higher resolution runs is A a 
twice as good as the standard resolution. 
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